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\q A measurement of the lifetime of the tau lepton has been made using a sample of 1671 Z° — > 
^ t + t~ decays collected by the SLD detector at SLC. The measurement benefits from the small 
and stable collision region at SLC and the precision pixel vertex detector of SLD. Three analysis 
techniques have been used: decay length, impact parameter, and impact parameter difference 
^ methods. The combined result is r T = 297 ± 9(stat.) ± 5(syst.) fs. 
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1 Introduction 



Measurements of tau lepton decays provide a unique test of lepton universality. Assuming the 
tau to have the same coupling to the W-boson as the muon and electron, the tau mass ra T , the 
tau lifetime r T , and the r — > ev e v T branching fraction B e are related as follows: 



where and are the mass and lifetime of the muon. Currently the precision of this test is 
limited by measurements of the lifetime and the electronic branching ratio of the tau. 

The Stanford Linear Collider (SLC) and the SLC Large Detector (SLD) provide an ex- 
cellent facility for measuring the tau lifetime. Decays of tau pairs produced at the Z° mass 
peak are highly boosted and collimated, and are relatively easy to distinguish from other final 
states. In addition, the SLC is characterized by a very small and stable luminous region, while 
the SLD with its CCD pixel vertex detector has the capability of measuring three-dimensional 
space points with high resolution close to the interaction point. The combination of these 
features substantially reduces the uncertainty in the determination of both the production and 
decay points of the tau, allowing for precise measurements with a relatively small sample of 
events. 

The results reported here are based on a sample of 1671 tau-pair candidates collected by 
SLD in 1992 and 1993. Three different techniques are used to determine the tau lifetime. In the 
first method, which uses tau-pair events in the 1 vs. 3 topology, the tau lifetime is extracted 
directly by measuring the decay length of the tau on the three-prong side of the event. The 
other two methods employ events in which both taus decay to a one-prong. In the impact 
parameter method, the tau lifetime is inferred from the distribution of impact parameters of 
the charged tracks, while in the impact parameter difference method, it is extracted from the 
correlation between the impact parameters and acoplanarity of the two tracks in the event. 
The decay length method gives a direct measurement of the tau lifetime, with relatively small 
backgrounds. The impact parameter and impact parameter difference techniques benefit from 
the large one-prong branching fraction of the tau, and are independent of the decay length 
method. 



A detailed description of the SLD detector can be found in Ref. [1]. A subset of detector 
elements relevant to the analyses reported here are described briefly below. These include the 
vertex detector (VXD) [2], the central (or barrel) drift chamber (CDC) [3], and the lead/liquid- 
argon calorimeter (LAC) [4]. 

The VXD consists of 480 charge-coupled devices (CCDs) surrounding a 1 mm thick beryl- 
lium beam pipe with an inner radius of 25 mm. Each CCD is an array of 375 X 578 square 
pixels 22 /mi on a side. The active material is a 20 fim thick epitaxial silicon layer on a 180 fim 
thick silicon substrate. The CCDs are mounted on 60 alumina boards 9.2 cm long, arranged 
in four concentric cylinders at radii ranging from 2.9 cm to 4.1 cm. The inner (outer) cylinder 
covers a range of polar angles defined by |cos#| < 0.85 (0.75). The CCDs are arranged so 
that at least two hits are possible over the full azimuth within the polar angle acceptance. On 
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average 2.3 CCDs are traversed by a track from the interaction point. The spatial resolution 
of the VXD is 6 //m transverse to the beam and 7 fim along the beam direction[6]. 

The barrel drift chamber in SLD is a cylinder 1.8 m long with an inner radius of 0.2 m and 
an outer radius of 1.0 m, filled with a CCvbased gas mixture. There are 640 drift cells arranged 
in ten superlayers covering radii from 24 cm to 96 cm, where superlayers with axial wires 
alternate with pairs of stereo layers at angles of ±43 mrad. Each cell in a superlayer has eight 
sense wires spaced radially by 5 mm. Field-shaping wires in each cell provide the desired drift 
fields. Each sense wire provides a measurement of the drift distance with a spatial resolution 
averaging 70 fim over the entire drift cell. Tracks are reconstructed with high efficiency at polar 
angles in the range | cos#| < 0.85. 

The LAC consists of an assembly of rectangular lead plates separated by 2.75 mm gaps and 
mounted in large insulated vessels filled with liquid argon. A barrel section covers polar angles 
in the range | cos#| < 0.82 and endcap sections complete the coverage down to | cos#| = 0.99. 
The LAC is segmented in depth into an electromagnetic (EM) section and a hadronic section 
(HAD). The EM section is made with 2 mm thick lead plates for a total depth of 21 radiation 
lengths and 0.84 hadronic interaction lengths, while the HAD section has 6 mm thick plates for 
an additional depth of 2.0 interaction lengths. The lead plates are segmented and connected so 
as to form projective towers with an azimuthal segmentation of 33 mrad in the EM section and 
66 mrad in the HAD section, and with comparable segmentation in polar angle. The energy 

resolution for electromagnetic showers has been measured to be about ^# = — , 15% 

to E y/E(G e y) 

The SLD event trigger requires any of several combinations of tracking and energy-flow 
information from the detector elements. A subset of these have a relatively high efficiency for 
tau-pair events, in particular a requirement of two back-to-back tracks, or a single track plus a 
minimum energy deposition in the calorimeter. 

2.1 Tracking Performance 

Charged tracks found in the CDC are linked to clusters of pixels in the VXD by extrapolating 
each track and selecting the best set of associated clusters [6]. A set of clusters may not be 
shared by multiple tracks. The track parameters are then recalculated, accounting for multiple 
scattering. In tau-pair events, at least one VXD cluster is linked to a well-measured CDC track 
in 99% of the cases. 

The measured impact parameter resolution transverse to the beam line for the combined 
tracking system is 11 /mi for infinite momentum tracks and 76 fim at 1 GeV/c. Along the 
beam direction the resolution is 38 fim for infinite momentum and 80 fim at 1 GeV/c. The 

momentum resolution of the tracking system in the 0.6 T field of the SLD solenoid is {^^j = 
0.01 2 + (0.0026 px) 2 , where px is the track transverse momentum in GeV/c. 

2.2 Interaction Point Determination 

The SLC collides bunches of electrons and positrons accelerated in the SLAC linac at a rate of 
120 Hz. After colliding at energies of 45.6 GeV each, the bunches are extracted and dumped. 
For the 1993 run, the spatial extent of the bunches at the interaction point was typically 
0.8 fim vertically, 2.6 fim horizontally, and 700 fim longitudinally [5]. The transverse 
position of the SLC collision region was stable, with variations of typically 5-10 fim over time 
periods measured in hours. 
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The spatial location of the interaction point is determined accurately in the transverse 
plane using tracks from hadronic Z° decays [6]. A fit to a single point in the transverse plane is 
made using tracks from 30 successive events, in an iterative procedure. At each iteration, tracks 
that significantly reduce the quality of the fit are removed. Typically the fit uses about 300 
tracks, and converges after a few iterations. The uncertainty in this determination is 7 ± 2 /mi 
for the 1993 data and 9 ± 2 /mi for the 1992 data. Non-Gaussian tails in the interaction point 
distribution may be represented conservatively by a second Gaussian with a standard deviation 
of 100 /mi for 0.25% (0.5%) of the 1993 (1992) events. 

3 Event Selection 

The events used in the present analyses were selected from a sample of 2.20 pb -1 collected at 
a center-of-mass energy of 91.2 GeV in 1992 and 1993. The selection of tau-pair candidates 
is based on the multiplicity, momentum and direction of tracks in the central drift chamber, 
and on properties of electromagnetic showers in the calorimeter. Tracks and energy clusters 
are required to meet certain criteria[7] to be used in the event selection. 

3.1 Selection Criteria 

Tau-pair candidates are required to have at least two but fewer than seven tracks. Each event is 
divided into hemispheres by the plane normal to the track with the highest momentum. Tracks 
in each hemisphere must fall within 15° of the net momentum vector in the hemisphere, and the 
jet invariant mass[8] in each hemisphere is required to be less than 2.3 GeV/c 2 . Furthermore, the 
jet axes in the two hemispheres must be back-to-back within 20°. These criteria discriminate 
strongly against background from multi-hadron final states. The polar angle m4 - ss of the 
missing momentum[8] in each event is required to satisfy |cos0 mj - ss | < 0.88 to discriminate 
against two-photon interactions and Bhabha events. Two-prong events are required to have a 
minimum acolinearity of 10 mrad. The scalar sum of the momenta of the two stiffest tracks 
in any event must be less than 65 GeV/c. These cuts primarily reject Bhabha events and 
muon-pair final states. 

The total visible energy[8] in an event is required to be at least 12% of the center-of-mass 
energy {Ecm) to reject two-photon interactions. To discriminate against Bhabha events, the 
total energy deposited in the electromagnetic section of the calorimeter is required to be less 
than 0.45 Ecm, an d the most energetic EM cluster must be less than 0.33 Ecm- In addition, 
the total calorimeter energy not included in identified jets[8] is required to be less than 5 GeV, 
and there must be fewer than six energy clusters not included in jets. 

These criteria resulted in a sample of 1671 tau-pair candidates selected from the 1992 and 
1993 data. Further requirements were imposed in each individual lifetime analysis. 

3.2 Selection Efficiency and Backgrounds 

The event selection efficiency and background contamination were estimated using Monte Carlo. 
The production of tau-pair events at the Z° resonance was simulated using the KORALZ 3.8[9] 
Monte Carlo generator. The same program was used to generate muon-pair events, while 
wide-angle Bhabha scattering, two-photon interactions, and Z° — > hadrons final states were 
produced using the generators described in Refs. [10], [11], and [12], respectively. All these 
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Monte Carlo data samples were subjected to the SLD detector simulation based on the GEANT 
3.15[13] program and to the above event selection. The SLD trigger was also simulated in the 
Monte Carlo. 

We determine a tau-pair event selection efficiency of 63.0 ± 0.3%, where the largest loss is 
due to the solid angle coverage of the CDC. The trigger efficiency is found to be 95.1 ±0.3%, and 
the result is an overall selection efficiency of 59.9±0.4%. Based on these estimates, the measured 
integrated luminosity, and the Z° — > t + t~ branching fraction, we expect a final tau sample of 
1714 ± 50 events, including backgrounds. This is consistent with the 1671 tau-pair candidate 
events selected in the data. From the fraction of Monte Carlo background events which survive 
the selection requirements, the purity of the tau sample is estimated to be 96.1 ± 0.4%. The 
accuracy of these estimates was checked by comparing data and Monte Carlo distributions for 
various quantities used in the event selection. 

3.3 Selection of Specific Topologies 

1-3 events: The decay length method (DL) described below uses only events with a 1-3 topol- 
ogy. Any pair of oppositely charged tracks which are consistent with originating from a photon 
conversion [6] are excluded. An event is then required to have exactly three tracks which form 
a common vertex in one hemisphere, and a single charged track in the opposite hemisphere. 
Each of the three tracks must have at least 25 hits in the CDC and at least one hit in the VXD, 
and the % 2 per degree of freedom of the track fit must be less than 5 for two of the tracks and 
less than 15 for the third. The % 2 probability for the vertex fit must be at least 0.02%. This 
results in a final sample of 257 events. 

1-1 events: For the two methods employing impact parameters, only events with a 1-1 topology 
are used (here photon conversions are also removed). For the direct impact parameter (IP) 
method, at least one track in each event must have: a momentum greater than 3 GeV/c and 
less than 40 GeV/c, at least 40 CDC hits and at least one hit in the VXD, a polar angle in the 
range | cos#| < 0.72, a % 2 per degree of freedom for the track fit less than 5.0, and a distance 
of closest approach to the interaction point along the beam direction less than 2.5 mm. For 
the impact parameter difference (IPD) method, both tracks in each event must satisfy these 
criteria. For both methods the tracks must have opposite charge, each event must have a two- 
prong invariant mass of at least 8 GeV/c 2 , the angle between the two tracks must be at least 
2.8 rad, and the missing momentum in the event must satisfy | cos m4 - ss | < 0.80. These criteria 
result in a sample of 1556 tracks (from 912 events) for the IP method, and 642 events for the 
IPD method. 



Process 


Backg 


;round fract 


ion (%) 




DL 


IP 


IPD 


Multihadron 


0.3 ±0.1 






Muon-pair 




0.9 ±0.1 


0.2 ±0.1 


Bhabha 




0.6 ±0.1 


0.1 ±0.1 


Two-photon 




1.0 ±0.3 


0.5 ±0.3 


Total 


0.3 ±0.1 


2.5 ±0.3 


0.8 ±0.3 



Table 1: Background fractions (in percent) in the final event samples 
used in the DL, IP, and IPD analyses. 

The main contamination for 1-3 events is from multihadron events, whereas for events in 

the 1-1 topology the primary backgrounds are Bhabha scattering, two-photon interactions, and 
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muon-pair final states. Table 1 summarizes the background fractions in the final event samples 
used in each of the three tau lifetime measurements. 




Three-prong Momentum (GeV/c) Three-prong Mass (GeV/c 2 ) 




7963M Opening Angle (deg) Vertex Probability 

Figure 1: Distributions from the 1-3 sample: (a) three-prong momentum, (b) three-prong invariant 
mass, (c) largest opening angle between a track and the three-prong momentum direction, and (d) 
three-prong vertex fit probability. In (d) the inset shows an expanded view of the lowest 0.4% in 
probability. 

A few track and event quantities are plotted in Fig. 1 and Fig. 2 for tau-pair events in the 
1-3 and 1-1 topologies, respectively. In all plots, the dots represent data and the histograms 
Monte Carlo. The three-prong momentum and invariant mass for 1-3 events are shown in 
Fig. 1(a) and 1(b), respectively. Figure 1(c) shows the largest opening angle between the 
three-prong momentum direction and one of the three tracks, while in Fig. 1(d) is plotted the 
probability of the three-prong vertex fit. The distributions in Fig. 2 for 1-1 events represent 
(a) the track momentum, (b) the track fit % 2 per degree of freedom, (c) the invariant mass 
of the two tracks and (d) their acolinearity. Good agreement is seen in all the distributions, 
indicating sufficient accuracy in the Monte Carlo simulation of the data. 
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Figure 2: Distributions from the 1-1 sample: (a) track momentum, (b) track fit x 2 /d.o.f., (c) two- 



prong invariant mass, and (d) two-prong acolinearity. 

4 Lifetime Measurement 

The three techniques used in this analysis are illustrated in Fig. 3. Figure 3(a) shows the princi- 
ple of the decay length method. Using tau-pair events with a 1-3 topology, selected as described 
above, a vertex fit in three dimensions is performed on the three-prong side to determine the 
decay point of the tau. Because of the relatively large beam size in the longitudinal direction, 
a two-dimensional decay length is calculated for each event using the precisely measured trans- 
verse beam position and the projection of the three-prong vertex in the transverse plane. The 
result is then translated into a three-dimensional decay length using the event thrust axis. 

Figure 3(b) illustrates the impact parameter method. In this technique, events in the 1-1 
data sample are used. For each event, either one or two measurements are extracted according 
to whether one or both of the two charged tracks in the event satisfy the selection requirements 
described in section 3.3. The impact parameter d of a daughter track is related to the transverse 
decay length l xy of the parent tau as follows (see Fig. 3(b)): 

d = l xy sin cf> = ! sin 6 sin </>, 
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where cf> is the angle between the track and the direction of the parent tau, and 6 the polar 
angle of the tau direction. 




Figure 3: Drawings illustrating the three measurement techniques used: (a) decay length, (b) impact 
parameter, and (c) impact parameter difference. The x and y axes lie in the transverse plane. 

The impact parameter difference technique [14] is illustrated in Fig. 3(c). Using events 
with a 1-1 topology, selected as described in section 3.3, this method exploits the correlation 
which exists between the difference of the impact parameters of the two tracks in the event and 
their acoplanarity. This is expressed by the following equation: 

< d + — gL > = < l xy > 6cf) =< / > sin# 6<f> } 

where d + and tf_ are the impact parameters of the two tracks, 6<f> is the acoplanarity of the 
tracks, and 6 is the polar angle of the tau direction. 

Directly or indirectly, all three methods involve the measurement of an average decay 
length < / >. The tau lifetime r T is related to < / > by: < / > = < ^7 > cr T , where < (Sj > 
is the average boost of the tau and is determined from Monte Carlo. In all three methods, the 
tau direction is taken to be the event thrust axis, which is determined on an event-by-event 
basis from charged tracks and isolated energy clusters in the calorimeter. Monte Carlo studies 
show this to be a very good approximation. 

In the decay length (DL) method, the result is derived from an unbinned maximum 
likelihood fit to the decay length distribution using an analytical function given below. In 
both the impact parameter (IP) and impact parameter difference (IPD) techniques, we make 
use of a binned maximum likelihood fit using Monte Carlo to parameterize the experimental 
quantities that are used in the two measurements: impact parameter for the IP method, and 
impact parameter difference and acoplanarity for the IPD method. 
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4.1 Decay Length Method 

The distribution of measured decay lengths for the 257 selected events is plotted in Fig. 4(a). 
The average decay length is extracted from a maximum likelihood fit using an exponential 
decay distribution function convoluted with a Gaussian resolution function: 

f f { x\ ( (x — U) 2 \ 
P(U, <r,-; Z , s ) = 7= / exp -- exp 1 dx, 

where the /; are the measured decay lengths, l is the parent decay length, and s is a scale 
factor on the calculated decay length errors <r 4 -. The Oi are determined on an event by event 
basis from a projection of the three-dimensional vertex error matrix along the event thrust axis, 
combined with a small contribution from the uncertainty in the transverse beam position. The 
mean value of the decay length error for these events is calculated to be 475 ± 9 /mi. From the 
fit, the average decay length is found to be 2.19 ± 0.14 mm, with a scale factor s = 1.31 ±0.11. 
The solid curve in Fig. 4(a) represents the fit function with these values, normalized to the 
data points. As a measure of the goodness of the fit, the % 2 per degree of freedom for the 
normalization is 0.8. 



600 




0.5 1.0 1.5 0.5 1.0 1.5 

Decay Length (cm) Decay Length (cm) 7963A4 

Figure 4: Three-prong decay length distribution for the final event samples in (a) the data and (b) 
the Monte Carlo. The solid curve in the two figures corresponds to the maximum likelihood fit 
described in the text. 

Figure 4(b) shows the decay length distribution for 2846 Monte Carlo tau-pair events 
generated with a mean lifetime of 305 fs and passing the same event selection criteria as the 
data. The fit yields an average decay length l = 2.331 ± 0.044 mm, corresponding to a mean 
lifetime of 305.6 ± 5.8 fs, in good agreement with the generated value. The average decay 
length error for these events is 460 /mi, and the fitted scale factor is s = 0.96 ± 0.05. Since 
the maximum likelihood fit to the Monte Carlo events returns the generated lifetime, no bias 
is attributed to the method and no associated systematic error is assigned. 

A Monte Carlo study using a large number of event samples showed that the fitted value 
of the scale factor s from the data is consistent with a statistical fluctuation from unity, as 



9 



the distribution of s is skewed toward larger values. As can be seen in Fig. 4, l is sensitive 
primarily to events with relatively large decay length, while s is determined mainly from decay 
lengths near zero, so that the lifetime result is fairly insensitive to the value of the scale factor. 
This was checked in the data by repeating the maximum likelihood fit with a fixed scale factor 
of unity. From the corresponding change in the fitted average decay length, a conservative 
systematic error of 0.8% in the lifetime is assigned. 

The effect of the track quality requirements and the cut on the vertex fit probability was 
studied in the data. From the changes in the result when these quantities were varied, we 
derive a systematic error of 1.0%. The sensitivity of the measurement to the event selection 
requirements was studied in the Monte Carlo, and the systematic uncertainty was found to 
be negligible. As seen in Table 1, the only significant background found to survive all the 
analysis cuts comes from multihadron final states at a level of about 0.3% of the 1-3 sample. 
It was determined that the average lifetime of such events is consistent with zero, and that 
including this background reduces the fitted decay length by 0.3%. Due to the uncertainty in 
the modeling of low-multiplicity hadronic events, a conservative systematic error equal to the 
full amount of this correction is assigned. 



Systematic effects 


Error (%) 


Decay length resolution 


0.8 


Track and vertex quality cuts 


1.0 


Backgrounds 


0.3 


Initial and final state radiation 


0.3 


Beam energy and energy spread 


0.3 


Total 


1.4 



Table 2: Systematic errors for decay length method. 



To study the effect of a possible misalignment between the VXD and CDC, the data were 
divided into four samples by azimuthal quadrants of the detector where decay vertices were 
found. Lifetimes determined separately for the four samples were the same within statistical 
errors, and the associated systematic uncertainty is estimated to be negligible. The effect of 
non-Gaussian tails in the distribution of beam positions has been estimated to be negligible in 
this measurement. The two-dimensional projection of the decay length using the thrust axis 
was checked by using instead the net momentum of the three-prong to define the tau direction, 
and results were found to be consistent. The uncertainty in the calculation of initial- and final- 
state radiation in the Monte Carlo was estimated to contribute a systematic error of 0.3% in 
the average boost of the taus. Finally, the effect of uncertainty in the SLC beam energy and 
energy spread was studied by varying these quantities in the Monte Carlo for tau events and 
observing the change in the average tau momentum. A systematic uncertainty in the lifetime 
of 0.3% is estimated from this source. 

The systematic effects discussed above are summarized in Table 2. From the Monte Carlo, 
we find for the average boost of the tau in these events: < f3~f >= 25.44 ± 0.01. Applying this 
factor, and including background correction and systematic errors, the decay length method 
yields a tau lifetime 

r T = 288 ± 18 ± 4 fs. 

This measurement was checked by an independent analysis in which the decay length was 
computed in three dimensions and considerable attention was paid to reducing the errors on 
the longitudinal position of the interaction point [15]. A consistent result was obtained. 
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4.2 Direct Impact Parameter Method 



As described in section 3.3, the event selection yielded a total of 1556 tracks for this lifetime 
determination. The distribution of impact parameters measured for these tracks is shown in 
Fig. 5. The impact parameter is assigned a positive (negative) value if the extrapolated track 
crosses the event thrust axis before (after) its point of closest approach to the interaction point. 
Negative impact parameters result from finite tracking errors and uncertainties in the beam 
position determination. 

To extract a lifetime from the impact parameter distribution, a binned maximum like- 
lihood fit to the data is performed. The fit function is represented by the impact parameter 
distribution from the Monte Carlo, corrected for background and normalized to the number of 
events in the data. The likelihood probability is expressed as follows: 

N 

In L = J2YMM0), 



i=i 



where N is the number of bins in the distribution used in the fit, Y{ is the number of entries 
in the i th data bin, and fi(9) is the normalized content of the i th bin in the Monte Carlo. The 
bin width chosen is 10 /mi, except that in the tails of the Monte Carlo distribution the bins are 
widened as required to include at least 10 events. The data are then binned the same as the 
Monte Carlo. 
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Figure 5: Impact parameter distribution for data (data points) and 
Monte Carlo (histogram). 

A single Monte Carlo sample was generated with an input lifetime of r = 305 fs. A 
weighting technique is then used to generate the impact parameter distributions corresponding 
to alternative lifetimes. A track originating from a tau decay with proper time t is taken to 
have been produced from a sample with different lifetime r by weighting its assigned probability 
with the following ratio: 

r exp(-t/r) 



weight 



t exp(-t/r ) 
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Weighted impact parameter distributions are formed with seven different lifetimes in the range 
0.7t o < t < 1.3r . These distributions are corrected for background by adding the normalized 
impact parameter distribution for background tracks passing all cuts. A likelihood probability 
is calculated for each distribution, and a fourth-order-polynomial fit is performed to the seven 
values as a function of lifetime. The best fit lifetime is taken to be the value corresponding 
to the maximum of this curve, with a statistical error assigned by taking the values where the 
likelihood has decreased by 0.5 from the maximum. The result is 

r T = 302 ± 12 fs. 

The impact parameter distribution corresponding to the best fit lifetime in the Monte Carlo, 
including the background correction, is shown as the solid histogram in Fig. 5. The agreement 
between data and Monte Carlo is very good, and the comparison between the two distributions 
gives a % 2 of 0.9 per degree of freedom. 

This analysis was checked for systematic bias using several different Monte Carlo samples 
of 1000 tracks each in place of the data. In every case the generated lifetime was reproduced 
within the statistical errors. A similar check was made with several groups of Monte Carlo 
events using the true generated tau direction in place of the event thrust axis, and results were 
consistent within statistics. The statistical error obtained from the likelihood fit was found to 
be consistent with the fluctuations in the results from the different Monte Carlo samples. 

The sensitivity of the fitting technique to bin size and Monte Carlo statistics was studied 
extensively. The analysis was repeated with nominal bin sizes a factor of two larger and smaller 
than the chosen 10 /tm, and the minimum number of events required in a Monte Carlo bin was 
also varied. The corresponding variations in the result lead to a systematic uncertainty of 0.3%. 
A study using several Monte Carlo samples of different sizes lead to an estimated systematic 
error of 0.9% due to finite Monte Carlo statistics. The range of impact parameters covered by 
the fit was varied, equivalent to a symmetric trimming of the sample from zero to 5%, both for 
the data and Monte Carlo. The results fluctuate within statistical errors, and no associated 
systematic error is assigned. 



Systematic effects Error (%) 

Binning 0.3 

Monte Carlo statistics 0.9 

Track quality cuts 1.0 

Backgrounds 0.1 

Beam position 0.4 

Tau branching ratios 0.2 

Initial and final state radiation 0.3 

Beam energy and energy spread 0.3 

Total 1.5 



Table 3: Systematic errors for the impact parameter method. 

All track quality cuts were varied in order to study their effect on the measurement. 
From the observed changes in the measured lifetime, we assign a conservative systematic error 
of 1.0%. The uncertainty in the background correction was checked by varying the background 
fraction in the Monte Carlo, and a corresponding systematic error of 0.1% is assigned. 

The possibility of detector misalignment was checked as described above for the decay 
length method, and no statistically significant effect was observed. The effect of non-Gaussian 
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tails in the beam position determination was studied in the Monte Carlo by introducing an 
additional smearing in the beam position in a fraction of events. The change in the measured 
lifetime was small and a conservative systematic error of 0.4% is assigned. We have also 
investigated the effect of uncertainties in the tau one-prong branching ratios used in the Monte 
Carlo. We estimate a systematic error of 0.2% from this source. Finally, as mentioned above 
for the decay length method, we assign a systematic error of 0.3% due to initial- and final- 
state radiation, and another 0.3% due to the uncertainty in the beam energy and beam-energy 
spread. 

A summary of the systematic errors is given in Table 3. Including systematic errors, the 
direct impact parameter method yields a tau lifetime 

r T = 302 ± 12 ± 5 fs. 
4.3 Impact Parameter Difference Method 

As described in section 3.3, the event selection yielded a sample of 642 events for this analysis. 
The impact parameter difference and acoplanarity distributions for these events are plotted in 
Fig. 6(a) and 6(b), respectively [16]. The data (represented by the dots) and the Monte Carlo 
(by the histograms) are seen to be in good agreement. In Fig. 7(a), the scatter plot of impact 
parameter difference versus acoplanarity for the data is shown. A clear correlation is evident, 
and the tau lifetime can be extracted from the slope of the distribution. This is even better 
illustrated in Fig. 7(b) where each data point represents the mean impact parameter difference 
< d + — d- > over a given bin in sin# 6cf>. 




-0.05 0.05 -0.2 0.2 

7963A6 d + - d_ (cm) sin6 8(|) (rad) 

Figure 6: (a) Impact parameter difference and (b) acoplanarity distributions for the final event 
sample in the IPD method. The dots represent the data and the histograms the Monte Carlo. 

A straightforward way to extract the tau lifetime is to fit the distribution in Fig. 7(b) 
to a straight line over a range in acoplanarity where the correlation is linear, together with a 
trimming procedure in order to remove from the fit tails that are not well modeled. In our case, 
we have chosen an alternative method based on a binned maximum likelihood method. This is 
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motivated by the fact that in any acoplanarity bin, the impact parameter difference distribution 
is asymmetric and non-Gaussian (it is exponentially distributed) and any truncation would 
result in a bias in the lifetime. 

The maximum likelihood technique used here is similar to that described in section 4.2 
for the single-impact parameter method; the only difference is that the data are described in 
two-dimensional bins in impact parameter difference and acoplanarity. Just as in the impact 
parameter method, a single Monte Carlo sample was generated, and several samples with 
different lifetimes were simulated using a weighting technique with a weight defined by: 



weight = — 



Tq exp 



(ti+t 2 ) 



exp 



(h +t 2 ) 

TO 



where r is mean lifetime in the original Monte Carlo sample, r the alternative lifetime, and t\ 
and t% the decay times of the two taus in the event. 
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Figure 7: Impact parameter difference vs. acoplanarity scatter plot for 1-1 events in the data. The 
points outside the parallelogram in (a) are excluded from the fit. The distribution in (b) represents 
a profile histogram of the scatter plot for all points used in the maximum likelihood fit. 

Before performing the fit, ten events lying in the tails of the impact parameter and 
acoplanarity distributions were removed (this corresponds to a 1.4% trim). This was done by 
defining a fit region represented by a parallelogram surrounding the core of the scatter plot in 
Fig. 7(a). In addition, in order to reproduce the distribution of impact parameter difference 
vs. acoplanarity as provided by Monte Carlo, the binning used in the maximum likelihood fit 
was chosen to be very fine in the middle of the scatter plot and progressively coarser towards 
the tails. In order to account for background, Monte Carlo events from the various sources of 
contamination were merged with the tau Monte Carlo sample used in the fit. The result of the 
fit for the tau lifetime using this technique is: 

r T = 298 ± 13 fs. 

The error is statistical only. A comparison between the scatter plot in Fig. 7(a) and its coun- 
terpart from Monte Carlo yields a % 2 per degree of freedom of 1.2. 

The largest source of systematic errors comes from the fitting procedure. We assign a 
systematic error of 1.6% due to the sensitivity of the measurement to the chosen bin size. This 
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is due mainly to the tails of the impact parameter difference and acoplanarity distributions 
where the finite statistics in the Monte Carlo do not allow a very fine binning. In the Monte 
Carlo sample, we require a minimum of 5 events in each square bin in impact parameter 
difference and acoplanarity in order to compute the likelihood probability for that bin. When 
this is not satisfied, entries from neighboring bins are progressively combined until the minimum 
number of events is reached. We have varied the number of entries required from f to f and 
observe no change in the measured lifetime. We have also studied the effect of the size of the 
fit region which we have varied over a wide range. This lead to assigning a systematic error 
of f .0%. Furthermore, we estimate a systematic error of 0.7% due to Monte Carlo statistics, 
determined from the spread in the measured lifetime as a function of the size of the analyzing 
Monte Carlo sample used in the maximum likelihood fit. 
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Figure 8: Impact parameter difference vs. 
wide-angle Bhabha events from the data. 
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Various analysis cuts were varied to study the effect of the event selection and the track 
quality criteria. We observe a 0.9% change in the measured lifetime when we vary the cuts 
that had the most effect at reducing the remaining background in the f-f sample, namely the 
cuts on the two-prong mass, the maximum track momentum, and cos0 m8SS . By varying the 
minimum track momentum and the minimum number of VXD and CDC hits on a track, we 
assign a systematic error of 0.5% due to track quality cuts. The sensitivity to the detector 
alignment was studied by dividing the data into bins in both the azimuth and polar angle; no 
significant effect was observed. 



Systematic effects 



Error (%) 



Binning 1.6 

Fit range f.O 

Monte Carlo statistics 0.7 

Event selection 0.9 

Track quality cuts 0.5 

Tau branching ratios 0.2 

Initial and final state radiation 0.3 

Beam energy and energy spread 0.3 



Total 



2.3 



Table 4: Systematic errors for the impact parameter difference method. 
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Similarly to the IP method, we assign a systematic uncertainty of 0.2% due to possible 
variations in the impact parameter distribution associated with each one-prong decay channel 
of the tau. Furthermore, we estimate a contribution of 0.3% from both initial- and final-state 
radiation and the uncertainty in the beam energy and the beam-energy spread. 

Other systematic effects were investigated. The analysis was performed on data samples 
from different run periods; the results were all consistent. An important systematic check of 
the impact parameter difference technique was done by running on event samples with known 
zero lifetimes, namely Bhabha scattering and muon-pair final states. This is illustrated in 
Fig. 8 which shows the impact parameter difference versus acoplanarity for (a) muon-pair and 
(b) wide-angle Bhabha events, both taken from the data. The two distributions are flat as 
expected, and the majority of the events are clustered at the origin[17]. As a result, this 
method is relatively insensitive to background. 

Replacing the data sample by several Monte Carlo data sets with the same number of 
events resulted in a spread in the measured lifetimes consistent with the statistical errors, while 
their mean value was very close to the input value in the Monte Carlo. This is a check that the 
method does not bear any systematic bias. This was further confirmed using several generator- 
level Monte Carlo samples with different lifetimes; the value of the tau lifetime with which the 
sample was generated was reproduced in each case. 

The systematic errors are summarized in Table 4. Adding all contributions in quadrature 
leads to an overall systematic error of 2.3%. Thus, the result from the impact parameter 
difference method is: 

r T = 298 ± 13 ±7 fs, 
where the first error is statistical and the second is systematic. 

5 Summary and Conclusions 

We have measured the tau lifetime using three different techniques. The results from the three 
methods are consistent with one another. The decay length method is independent of the 
other two techniques since it uses a completely separate set of events. The impact parameter 
and impact parameter difference methods share about 70% of the events. However, because 
they make use of different information, they are not totally correlated. We have determined 
the correlation between the two techniques using Monte Carlo and following the procedure 
described in Ref. [18]. For the events that are common to the IP and IPD methods, we obtain 
a correlation factor of 53%. This results in an overall correlation of 37% between the two 
measurements. In this evaluation, common systematic effects are also taken into account. 

A combined result of r T = 300 ± ll(stat.) ± 5(syst.) fs is derived from the IP and IPD 
measurements. With the inclusion of the DL measurement, we obtain: 

r T = 297 ± 9(stat.) ± 5(syst.) fs. 

This result is in agreement with the present world average value [19]. Our measurement is 
currently limited by the relatively small size of our data sample; we expect a significant im- 
provement in the future as we anticipate a substantial increase in statistics, together with an 
associated decrease in the systematic error. A precise measurement of the tau lifetime at SLD, 
performed in a different environment relative to other experiments because of the small beam 
size and three-dimensional high-resolution vertexing, will represent an important contribution. 
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